^gjR-OUD-WOO 


/M 7 / 

J 6'£ 

10 f 


NASA Contractor Report 189645 


ADAPTIVE NONLINEAR POLYNOMIAL NEURAL 
NETWORKS FOR CONTROL OF BOUNDARY 
LAYER/STRUCTURAL INTERACTION 


B. Eugene Parker, Jr. 

Richard L Cellucci 
Dean W. Abbott 
Roger L. Barron 
Paul R. Jordan, III 

BARRON ASSOCIATES, INC 
Route One, Box 159 
Stanardsville, Virginia 22973-9511 

and 


00 








IT 

m 


'O 

rsj 

fU 

in 

1 

r-» 

m 


U 

O 

O' 

c 

r\j 

2 

3 

O 


r- 

m 

O 


H. Vincent Poor 

Department of Electrical Engineering 

Princeton University 

Princeton, New Jersey 08544-5263 


Contract NAS1-19271 
December 1993 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23681-0001 


V c 

Q C '~ 

< LL H 

O O' 

2 2 a 

JOQh 

uj < a ^ 

> or ca k * a 

3 U CT 

h- UJ UL < DO 

a 2 o a: < 

< uj 

0 J JhlO 

< < O 2 r-l 

H Y H 

X ID 

^ O ? J +- 1 

^ 2 a < * na 

^ V u a: c *- 

sO -J O ^ u 

ch o or -5 o 
foao u ui 

rH 11 D sO W 

1 QCT (Y H < 

or < oo k 

0 uj ^ co •> c 

1 2 a: \ v c 

< *-* a a: u w 

^ J3 u c w 

< 2: k > a -a 

Z C Ui < D CD 

^ 2 z -j esc: ^ 



Foreword 

This Phase I Small Business InnwaBon Research (SBIR) final report has been 
prepared by Barron Associates, Inc., Stanardsville, Virginia, to document the work 
performed from January 16 to August 15, 1991, under Contract NAS1-19271 with the 
NASA Langley Research Center, Structural Acoustics Branch. 

The authors express their gratitude to Drs. Richard J. Silcox, Technical Monitor, 
and Lucio Maestrello for their assistance in the formulation of the models utilized in 
this report and for their encouragement and support. 

Opinions expressed in this report are those of the authors and Barron Associates, 
Inc., who are solely responsible for its content. 


Abstract 


The acoustic pressures developed in a boundary layer can interact with an aircraft 
panel to induce significant vibration in the panel. Such vibration is undesirable due to 
the aerodynamic drag and structure-borne cabin noise that result. The overall 
objective of this work is to develop effective and practical feedback control strategies 
for actively reducing this flow-induced structural vibration. 

This report describes the results of initial evaluations using polynomial neural 
network-based feedback control to reduce flow-induced vibration in aircraft panels due 
to turbulent boundary layer/structural interaction. Computer simulations are used to 
develop and analyze feedback control strategies to reduce vibration in a beam as a 
first step. The key differences between this work and that ongoing elsewhere in the 
active control of structural vibration are, firstly, that turbulent and transitional 
boundary layers represent broadband excitation and thus present a more complex 
stochastic control scenario than that of narrowband (e.g., laminar boundary layer) 
excitation; and, secondly, that the proposed controller structures are adaptive 
nonlinear infinite impulse (HR) response polynomial neural networks, as opposed to 
the traditional adaptive linear finite impulse (FIR) response filters used in most 
studies to date. 

The controllers implemented in this study achieved vibration attenuations of 30 to 
60 dB, depending on the type of boundary layer established by laminar, turbulent, and 
intermittent laminar-to-turbulent transitional flows. Application of multi-input, 
multi-output, adaptive, nonlinear feedback control of vibration in aircraft panels 
based on polynomial neural networks appears to be feasible today. Plans are outlined 
for Phase II of this study, which will include extending the theoretical investigation 
conducted in Phase I, and verifying the results in a series of laboratory experiments 
involving both beam and plate models. 
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Introduction 


The acoustic pressures developed in a boundary layer can interact with an aircraft panel 
to induce significant vibration in the panel. Such vibration is undesirable due to the 
aerodynamic drag and structure-borne cabin noise that result. The overall objective of 
the proposed work is to develop effective and practical feedback control strategies for 
actively reducing this flow-induced vibration of aircraft panels. 

This report describes the results of initial evaluations using neural network-based 
feedback control strategies to reduce flow-induced vibration due to turbulent boundary 
layer/structural interaction. Computer simulations are used to develop and analyze 
feedback control strategies to reduce the vibration in a beam as a first step. There are 
two key differences between this work and that ongoing elsewhere in the active control 
of structural vibration. Firstly, that turbulent and transitional boundary layers repre- 
sent broadband excitation and thus present a more complex stochastic control scenario 
than that of narrowband, deterministic excitation, such as from propeller-induced noise. 
Secondly, that the proposed controller structures are adaptive nonlinear infinite im- 
pulse response (HR) polynomial neural networks, as opposed to the traditional 
adaptive linear finite impulse response (FIR) filters used in most studies to date. 

Simulation involved two basic modeling tasks: modeling the beam response and 
modeling the excitation field. Both of these issues are discussed herein. With regard to 
the panel response, the basic assumption invoked is that this response can be described 
as linear, as is conventionally done in studies of panel vibration (see, e.g., [25]). Such 
a model can be developed from the basic principles of plate theory [5, 8], and this is 
the approach that is adopted here. ... 

Due to the complexity of the flow processes of interest, useful models for the excita- 
tion field are necessarily stochastic. For self-preserving boundary layers, such as those 
occurring in purely turbulent flow, the excitation field can be assumed to be a Gaussian 
random field (e.g., [3]). Under this assumption, it is sufficient to model the excitation in 
terms of its second-order statistical properties, specified in terms of correlation or spec- 
tral properties of the excitation. Experimentally-derived models for the second-order 
structure of the turbulent boundary layer (as it interacts with the panel structure) are 
reported in [17-19, 24]. Certain of these models will be adopted here as the bases for 
the development of linear stochastic control stategies for flow-induced panel vibration 
due to turbulent boundary layers. 

In addition to studying the control of vibration in response to a turbulent boundary 
layer, this initial study also focuses on the control of vibration induced by laminar- 
to-turbulent transitional boundary layers. Transitions from laminar to turbulent flows 
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result in non-equilibrium boundary layers, which cannot be modeled accurately as 
Gaussian random fields. In particular, the behavior of the pressure at a point in a 
laminar-to-transition zone region resembles an intermittently gated random signal [5, 
p. 579], Even though the random signals themselves are Gaussian (i.e., they are the 
pressures of purely turbulent flow, or of purely laminar flow), the resultant of the 
intermittent gating produces a non-Gaussian signal. Modeling of this type of excitation 
for the purposes of vibration control can take one of two forms, depending on the control 
strategy to be applied. If a linear control structure is to be used, then it is sufficient to 
model the excitation based on its second-order statistical properties, specified in terms 
of the correlation or spectral properties of the excitation. Second-order models for the 
pressure field due to laminar-to-turbu!ent transition zones are described, for example, in 
[9] and [5, Eq. 8-71]. These models, perhaps with suitable modifications, can be used as 
the bases for developing linear stochastic control stategies for transitional-flow-induced 
panel vibration. 

Optimum linear feedback control based on second -order modeling is globally opti- 
mum for the control of Gaussian vibration fields. However, the non-Gaussian excitation 
found in transitional boundary layers will produce a non-Gaussian vibration field in 
the panel to be controlled, and the corresponding optimum stochastic control strategies 
will therefore be nonlinear. Effective and adaptive approximations to such optimum 
control strategies are also likely to be nonlinear, and as such can be implemented using 
the general structure of nonlinear polynomial networks (e.g., [4]). Since these transition 
zone regions are potentially a major source of vibration in aircraft panels, it is essential 
that their full statistical behavior (i.e., not only their second-order characteristics) be ap- 
propriately modeled in the process of developing vibration control strategies. Even if 
linear control strategies based on second-order models were to emerge as being nearly 
as effective as fully optimal nonlinear solutions (an eventuality that is not clearly likely), 
their relative effectiveness could only be assessed through consideration of higher-order 
statistical characteristics of the flow. 

This report is organized as follows. Section 2 describes the fundamental stochastic 
partial differential equations (PDEs) that model the response in the linear regime of 
a flat plate immersed in a fluid, when subjected to a boundary pressure field. Sec- 
tion 3 discusses the acoustic damping arising in such a model due to the interaction 
of the fluid with the plate. Two basic models are discussed in this context, corre- 
sponding to the cases of subsonic and supersonic fluid motion. Section 4 describes 
two experimentally-determined stochastic models, again corresponding to the subsonic 
and supersonic regimes, for a turbulent boundary convected over the panel. Section 5 
discusses the decomposition of the panel motion into its modal behavior, and Section 
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6 discusses the modeling of sensor and actuator configurations. In Section 7 , the over- 
all model is described in some detail for the lower-dimensional problem of controlling 
a vibrating beam, and in Section 8 a hierarchical succession of benchmark examples 
is proposed for computer simulation of the boundary layer/structural interaction. In 
Section 9, experimental results are given demonstrating the effectiveness of feedback 
control in active vibration suppression for a beam excited by laminar and turbulent 
boundary layers, as well as laminar-to-turbulent transitions of the boundary layer. Fi- 
nally, in Section 10, future directions for continuance of this work in Phase II are briefly 
outlined. 


2 A Fundamental Model for Panel Response 

We wish to model the motion of a tHin, flat rectangular panel occupying the region 
7Z = {0<;X ^ a i 0 < y < 6, z = 0}, over which a turbulent fluid is flowing at an 
average velocity v in the positive x direction. It is assumed that the plate is of uniform 
stiffness, and is clamped at the edges. 

The motion of such a plate is described by the stochastic PDE [8] 

|l)V 4 + j w = / + d, (2.1) 

where u>,/, and d are fields (in this case, functions of x,y, and time t) representing 
the plate displacement, the external pressure due to the turbulent boundary layer, and 
damping, respectively, all taken as being positive in the negative z direction. Here, V 4 
is the biharmonic operator 

V ._(®L + 2 _*l_ + » 1). ( 22 ) 

Ui< + d*VV + dy'j ' ' ’ 

D is the plate rigidity; and p\ is the mass density (per unit area) of the panel. The 
plate Eq. 2.1 is subject to the boundary conditions imposed by the clamped edges: 

iu(0, y,J) = = u;(x, 0, t) = w(x,b,t) = 0. (2.3) 

It is assumed that the external pressure field is due exclusively to a turbulent bound- 
ary layer on the z > 0 side of the plate. In this case, the forcing field / can be modeled 
as a homogeneous Gaussian random field whose second-order statistics will be de- 
scribed below. 
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The damping field d(x,y,t) can include both structural and acoustic damping 
terms. Structural damping is modeled by including a term = 7 j 

> 0 is a coefficient of structural damping [19]. Acoustic damping is due to the 
feedback effects of an acoustic velocity potential, <p(x,y,z,t ), launched into the fluid 
by the plate motion [24]. The evolution of this potential is governed by the wave 
equation 

c ! v ’*=(! + 4 ) * -’ >o - (2 ' 4) 

where V 2 = is the Laplacian operator, and c is the speed of sound 

in the fluid. This potential is driven by the plate motion via the boundary condition: 

d <9 

*,<)!.« o = -Q t w(x ,y,t). (2.5) 

The resultant damping field coupled into the plate is then given by 

d*{x,y,t ) - p 2 ^ <p{x,y,0,t), (2.6) 

where pi is the mass density of the fluid. 


3 Damping Models 

We will consider two specific damping models, corresponding to subsonic and super- 
sonic flow, respectively, within the general model described above. 

3.1 Subsonic Flow 

For subsonic flow (v < c), the viscous damping can be treated as an additional term of 
the form d a = -~f a §- ( w with 7 a > 0 [19]; thus the equation of interest is a second-order 
(in time) stochastic PDE: 

( CV « + ,,|l + 4) a ' = / ' (3-D 

^ ^ jhis type of damping is know as proportionate damping; and, 

although it is somewhat ad hoc, it has the virtues of parsimonious description and 
intuitive simplicity. Unfortunately, this model is too simple to accurately describe the 
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damping observed in physical plates. A related damping model that is more consistent 
with observed behavior introduces proportionate damping terms in each mode of the 
plate, allowing 7 to vary with mode number [19]. 

In particular, with proportionate damping, we note that the behavior of a given 
mode (say, the m th ) in the eigenfuntion decomposition of the panel motion will have 
the dynamics of a simple mass-spring system, described by the following scalar second- 
order ordinary linear differential equation: 

( p ' 7 £ + '‘It + DK ) Wm(t) = (3 ' 2) 

The dynamics of this system are those of a damped harmonic oscillator with damping 
ratio £ = 7/2^' m pi and damped resonant frequency u-y = u: m n/ 1 — £' 2 , where the 
undamped natural frequency w’ m equals \fD~^nJpi- A closer fit to actual panel be- 
havior can be achieved by allowing 7 to vary with m. That is, we consider temporal 
modal displacements w m (t) satisfying the differential equation 

+ BA -) “■■■>(<) = /"( () - < 3 - 3 ) 
This does not translate back into the simple global proportionate damping model of 
Eq. 3.1 when the modes are superimposed^ Rather it provides a richer damping model 
while preserving the linearity and anafytical tractability of the eigen-decomposed plate 
model. 

Of course, the individual modal dynamics are still those of mass-spring systems 
with the only modification being that the damping ratios become £ m = 7 m /2u; m p 1 . 
As with a global damping based upon 7, the individual modal damping coefficients 
7 m must be chosen empirically. For example, one choice suggested in [19] for a specific 
class of plates is 7™ oe Typically, the range of damping ratios £ that one would 

expect in panels are on the order of 0.001 to 0.01. 

3.2 Supersonic Flow 

For supersonic flow (u > c), the modeling of the acoustic damping is more complex. 
In this case, the relationship between the motion of the plate and the acoustic potential 
at the surface of the plate is described by [25]: 

^(x,y,0,uj) = jf G{x - x\y -y’,uj) + v-^J w{x’,y' ,u)dx’dy', 

(3.4) 
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where <j> and w are the Fourier transforms of Q and w, respectively, as functions of 
their time variable t; and where G is a Green's function for the system of Eqs. 2.4-2.6, 


given by: 


G(x,y,u>) = 


0 otherwise 


(3.5) 


where k = u :/c, and M is the Mach number [M = v/c). " 

The integral in Eq. 3.4 has a strong singularity of the Cauchy type along the inter- 
section of the Mach cone with the z = 0 plane. Thus, the evaluation of this transfer 
characteristic requires some care. 


4 Stochastic Models for the Turbulent Bound- 
ary Layer 

The forcing field / in Eq. 2.1 is the pressure exerted on the plate by the turbulent 
boundary layer. This field can be modeled as being a homogeneous Gaussian field [15], 
and thus its specification requires only the determination of the second-order correlation 
function of the field; viz.: 

T\(Z,ti,t) = {f{x,y,t)f*(x + £,y + T),t -1- r)). (4.1) 


We consider two different models for this correlation structure, corresponding to 
experimentally-derived models for the subsonic and supersonic cases. 

In the supersonic case, we consider the model [25]: 


<1 (OAK 1 


(4.2) 


where the As, I\ s and os are constants, 8 is the boundary layer thickness, U c is the 
convective velocity (assumed to be less than v), and (p^) is the mean-square intensity 
of the forcing field. Experimentally determined values for the A n and K n parameters 
are given by Ai = 4.4 x 10 2 ,A^ — 7.5 x 10 2 ,A 3 = —9.3 x 10 2 , A\ = 
-2.5 x 10~ 2 , I<i = 5.78 x 10“ 2 , K 2 = 2.43 x 10 ~\I< 3 = 1.12, and K< = 11.57. 
Specific values of the remaining parameters for a benchmark example will be given 
below. 
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An empirically-determined model for the pressure field in the subsonic case is 
described by the covariance structure [17]: 


riOe^r)^ 2 )^ 


lA'K'c-W** 


„tU( A ») 2 + (i /mu-UcTy + r,*} 


( 4 . 3 ) 


where 0 , the A's and I\ f s are constants; and where 8* denotes the boundary layer 
displacement thickness. Note that this is only one of several similar models that have 
been proposed for the pressure field of a subsonic turbulent boundary layer [17, 19], 
and it is chosen primarily because of its functional similarity to Eq. 4.2 in the two key 
variables (viz. £ and r). Experimentally determined values of A' n and I\ l n parameters 
are given by A\ = 0 . 24 , A 2 = 1 . 08 , A 3 = 1 . 80 , K[ = 0 . 47 , A 2 = 3 . 0 , and A 3 = 
14.0. As for the supersonic model, benchmark values of the remaining parameters will 
be specified below. 

As noted in the Section 1 , this Gaussian model for the pressure induced by the 
boundary layer is not suitable for transitional boundary layers. Later in this report, 
we implement a model of the transitional boundary layers, representing non-Gaussian 
excitation. To do so, we consider a model that treats the pressure field under all condi- 
tions as a mixture of those induced by purely laminar and purely turbulent flows, each 
of which can be modeled as a Gaussian random field. (Actually, the laminar flow can 
be modeled as being sinusoidal without significantly sacrificing accuracy.) The mixing 
of these two flows can be modeled by an independent, binary gating field, that switches 
intermittently between the two equilibrium boundary layers. This binary field can, in 
turn, be modeled as being Markovian in space and in time (i.e., it can be modeled as 
a dynamic Gibbs field), with transistion statistics that can be parameterized to yield 
either of the pure states or any degree of transition. This then, will add parameters 
to the model in addition to those discussed above. It should be noted that such a 
strategy could also be used to develop more refined second-order models, analogous 
to those of Eqs. 4.2 and 4.3, for the transition-zone statistics; although a full-order sta- 
tistical model is of primary interest in the present context of developing algorithms for 
vibration control. 


5 Evolution in Terms of Spatial Modes 

In order to create a more parsimonious desciption of the evolution described by Eq. 
2 . 1 , it is useful to represent all spatial characteristics in terms of the modes of the plate 
via a set of orthonormal eigenfunctions of the biharmonic operator, V 4 . In particular. 
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on denoting by {A*}^ and the eigenvalues and corresponding orthonormal 

eigenfunctions of V 4 (i.e., solutions to V 4 ^ = \ip), we can then write Eq. 2.1 as the 
vector differential equation: . 

if l dP + ^ v dt + DAlJ w(t) = /(<) + a(<), (5.1) 

where, A is the diagonal matrix with diagonal elements Aj, Aj, . . . ; and for each h ~ 
1, 2, .... the k th components of w(t), /(f), and c(t) are given by 


w k (t) = J w(x,y,t)ip k (x,y)dxdy, 

(5.2) 

fk(t)= f j f{x,y,t)xl> k {x,y)dxdy, 
Jo Jo 

(5.3) 

a*(0 = / / d a (x,y,t)ip k (x,y)dxd.y. 
Jo Jo 

(5.4) 


(Of course a k (t) is simply — laj- t c k{i) ' n the subsonic case.) The displacement field at 
each time t is then given by the expansion 

OO 

w{x,y y t) = v>k(t)i>k(x, y), (5.5) 

k= 1 

Further discussion of this expansion, including remarks on the structure of the input 
and damping terms, is included below. 

6 Sensor and Actuator Modeling 

In order to study the control of the field w we can consider three models of interest. In 
the first of these models, we assume that the field can be sensed throughout the region 
It, and that we can apply controls at all points in It. That is, we assume measurements 
of the form -t..,. - 


r(x,y,t) = w(x,y,t) + n(x,y,t), (x ,y)e1t, t > 0, (6.1) 

where n represents a white Gaussian measurement noise field; and a control input to 
Eq. 2.1 of the form 

c(x,y,t)=[ f f h(x,y,t\x ,y' ,t')r{x' ,y\t')dx'dydt' , (6.2) 

Jo Jo Jo 


8 


where h is the transfer function of the controller. Note that, in assuming a linear 
control law, there is no loss in generality with a quadratic loss performance criterion 
(e.g., minimizing squared error) since the field to be controlled is Gaussian and the 
measurements are linear. 

This model can be decomposed into modal components. In particular, we can 
consider the equivalent model: 


r(t) = w(t) + n(t), (6.3) 

where n represents a vector of independent and identically distributed Gaussian white- 
noise processes. Within this representation, the control signal becomes 

e(t) = £H{t,t')r{t’)dt', (6.4) 

where the matrix function H has k — I th element 

= [fff h(x,y,t;x\y',t')il} k (x\ y')xpi{x, y)dxdy'dxdy. 

J 0 Jo Jo Jo 

(6.5) 

This first model is primarily of theoretical interest since it is not practical to apply 
sensing and actuating materials to the entire panel. However, consideration of this 
model can provide useful insight into the controllability properties of the overall system 
describing the plate motion. 

As a second, more realistic measurement model, we can consider the sensing/actuation 
of the field by point sensors/actuators placed on a lattice of points, i = 

1, . . . ,n,j = l,...,m} in 71. In the simplest case of this model, we can assume 
measurements of the form 

r(x»yj,t) = w(xi,yj,t) + i = 1, . ... ,n, ; = (6.6) 

where i = l,...,n, j = 1, m} is a family of independent and 

identically distributed white-noise processes. This model can be modified to include 
the dynamics of the sensors by replacing the terms w(xi,ijj,t) with filtered versions 
w'(xi,yj,t) given by 

~T w '{ x i,yj,t) = -fsxv'{xi,y h t) + w{xi,yj,t), (6.7) 
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where f a > 0 determines the bandwidth of the sensors. (Obviously, this model assumes 
identical, first-order sensors.) In terms of the modal response, the model of Eq. 6.6 can 
be written as 

OO 

r(x„y,,f) = w k (t)ipk(xuyj) + i = 1,. •• ,n,i = I, . . . ,m; (6.8) 

k=o 

or to incorporate sensor dynamics, we can replace the w jts with w' k s generated by 

~ -f* w k ( 0 + MO- (6.9) 

Note that, within this model, we have multiple independent measurements of each 
component of r(f), but that these measurements are not independent from component 
to component Of course, as the lattice becomes denser, these various measurements can 
be combined to provide approximately independent measurements of the components 
by approximating the continuous model of Eq. 6.3. 

In the point-actuator model, the control signal will take the form 

n m 

c(x,y,t) = j]53c(x, -,%,*)£(* ~ *»)% ~ V})'- (x,y) € TZ.t > 0, (6.10) 

«'=i j=i 

where S denotes the Dirac delta function, and where the coefficients c(x„ yj,t) are 
obtained via linear transformation of the measurements for times up to time t\ that is, 

c(x, ,!/_,,<) = YY / ^ x i',yj'i t')dt' , (6.11) 

where the coefficent function h determines the controller. In terms of modal response, 
the vector control signal analogous to that of Eq. 6.4 has k th component in this case 
given by 

71 771 

<M) = Y Y c ( Zi ’ !/>> 0M*.> M (6.i2) 

;=1 j = 1 

Analogously to Eq. 6.7, actuator dynamics can be included in this model by replacing 
the signals c(x,-, yj, t) with filtered versions c'(x,, i/j, t) 

= -f a c'{x u yj,t) + c{xi,yj,t), (6.13) 
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where f a > 0 controls the actuator bandwidth. Of course, both Eqs. 67 and 6.13 can 
be incorporated into the plant when designing and evaluating specific control laws. As 
a final comment on this model, we note that this model is easily modified to allow for 
sensing and actuating on different lattices. 

The third model of interest is that in which the sensors and actuators are placed on 
lattice points as above, but they are modeled as having significant dimensions. In this 
case we can think of an observat ion for each lattice point that integrates the displacement 
field over the subset of the panel occupied by the sensor at that lattice point. That is, 
on denoting the region occupied by the i — j th sensor by 7 Zij, we can integrate the 
continuous model of (6.1) to yield the discrete model 

i = I, . . . ,n,j = 1, m. (6.14) 

where Rij(t), and -V i ; (f) denote the integrals of the fields r(x, y. t), w(x, y, <), 

and n(x. y, t ). respectively, from Eq. 6.1, over the region 7 Zij. Note that we retain the 
notation of Eq. 6.6 for the noise, since the noise elements A 'ij(t) obtained by integrat- 
ing n(x, y.t) over the regions 7£,j form a matrix of independent white noises as was 
the noise model in Eq. 6.6. These noises will be identically distributed if all sensors 
have the same area. (It is assumed, of course, that the sensors do not overlap one 
another.) 

The i — j th actuator can be assumed to apply uniform pressure over the region 
7 Zij. Thus, the control signal in this model takes the form 

71 T71 

c(x.y.t) = ^,52Cij(t)in i Jx,y), (x, y) e 71, t > 0. (6.15) 

»=i j=i 

where In, denotes the indicator function of 7 Z,j] and where the functions Cij(t) are 
obtained by linear transformations of the measurements Eq. 6.14 in a manner analogous 
to Eq. 6.11. 

In modal form, the sensor model Eq. 6.14 is 

OC 

R iA t )= Y^ w k{t)$k(i,j) + Nij(t), i = (6.16), 

with = Jn, j ‘ l Pk[x,y)dxdy, and the k th comp>onent of the modal control 

signal becomes 

^) = EEW(‘-i)- (6.17) 

»=ij=i 
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As in the point sensor/actuator model, the sensor and actuator dynamics can be 
introduced by replacing the quantities Wjj(f), Cij(t), and Wk(t ) with appropriately 
filtered versions. Also, as above, the model is easily modified to provide for non- 
colocated sensors and actuators. 

For the Phase I analysis of the vibration control problem, it will be assumed that 
sensing and actuation are performed in a continuously distributed fashion. The modal 
controllers will assume complete knowledge of modal displacements (i.e., ideal sensing) 
and the ability to effect unimodal actuation (i.e., ideal actuation). Consideration of this 
model can provide useful insight into the controllability properties of the overall system 
describing the plate motion. This is appropriate for the Phase I study of the viability 
of active control. 

We end this section with a few comments. First, note that any of the three sensing 
models can be combined with any of the three actuation models by a straightforward 
modification of the equation describing the transfer from measurements to controls. 
Secondly, it should be kept in mind that in a practical system, these transfer character- 
istics will take on a constrained form. (For example, their time behavior will likely be 
constrained to be that described by time-invariant, lump-parameter systems.) Thirdly, 
in the event that the forcing and damping functions can be modeled as autoregressive, 
moving-average (ARM A) processes, the optimum stochastic controller will take the 
form of an optimum state estimator (e.g., a Kalman filter), followed by an optimum 
deterministic controller. Finally, we recall that the assumption of linearity in the control 
law will not be reasonable for the generalizations of this problem in which the forcing 
function is a transitional field [5, 11]. This is because such fields are not Gaussian, and 
therefore do not yield linear structures as optimal controllers even under quadratic loss. 
Thus, in the case of transitional boundary layers, nonlinear mappings (such as polyno- 
mial networks) from the sensors to the actuators will replace the linear mappings of 
Eqs. 6.2, 6.11, etc. 

7 Consideration of a Vibrating Beam 

One difficulty with the panel model described in the preceding sections is that it is 
quite complex spatially. In order to reduce this complexity, it is useful to consider first 
a y = 0 "slice" of the problem. This simplification will still allow us to investigate 
many of the basic phenomena of interest, provided that the slice is oriented in the 
direction of fluid flow (i.e., in the x direction). Thus, we will consider the control of 
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the analogous stochastic PDE in two spatial dimensions (i and z): 

D— + pi-^ + 1s-Q t ) w = / + da + c ’ 0 < x < a ’ 1 - ( 7 . 1 ) 

where the displacement w, forcing field /, acoustic damping d a , and feedback control 
c are functions of x and t. The boundary conditions are 

w(0,t) = w(a,t) = 0. (7.2) 


In the supersonic case, we will consider the acoustic damping to result from an 
acoustic potential <t>(x, y, t) evolving in the x — z plane via the two-dimensional version 
of Eqs 2.4-2.6: 



The spatial eigenvalues and orthonormal eigenfuctions of Eq. 7.1 are given by (e.g.. 



where e* = ’ anc * ft, ft,--., are the roots of the equation cos /3 cosh (3 = 

1. These roots are given approximately by ft=4.730,ft=7.853, and ft=7r(2fc + 
l)/2, for k > 2. 

Knowledge of the eigenstructure of the clamped beam allows us to model the forcing 
function in the modal decomposition of Eq. 7.1. In particular, the forcing function can 
be written as 

f(x,t) = ( 7 - 8 ). 

Jt=i 
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with 


fk{t)= [ f{x,t)rj; k {x)dx, Ic = 1,2,.... 
Jo 


(7.9) 


Since / is a Gaussian field, the vector process /(f) will also be Gaussian, Thus, we 
can model its statistics by determining its autocorrelation matrix r(r), whose k — I th 
element is given by: 

r w ( r ) = (fk(t)fi{t + r)) = I i' k {x) [ Mx')r{x-x' t 0,7)dxdx'. (7.10) 

J 0 JO- 

where T is derived from Eq. 4.2 or 4.3, depending on whether we are considering the 
subsonic or supersonic case. In either case, the computation in Eq. 7.10 is simplified 
by first transforming to the temporal frequency domain. After Fourier transformation 
(in the r variable), both Eqs. 4.2 and 4.3 reduce to the form: 

f x (C,0,u ; )=>(u;)e- j|€l+ ^ /t/ 'r ' (7.11) 

where F is a weighted sum of exponentials; and where g = jjj in the subsonic case, 
and g = in the supersonic case. Thus, in the temporal frequency domain the 
elements of T are given by: 

fjuM ~ F{^) J° Mx) e ~ A ' X J Q M x ’)e A ‘*'dx' + e Ax j rJji(x)€~ Ax 'dx' dx, 

(7.12) 

with A = g + iu)/U c . By inserting Eq. 7.7 into Eq. 7.12, the forcing statistics can be 
determined in closed form. In particular, we have 

(7.13) 


f*,;(-') = F(w)(Ajt ) /(w) + A ( * tJt (u;) + d, fc,/(w) + — ^(w)], 


where 


A*.,f(<^) = 
and 


[e~ Aa (Av\ 2 \0) - 0 ; (3, (O)) - {A^\a) - tf°( a))] WEV) + 4 V)) 
’ (A 4 - A fc )(A 4 - A/) 

(7.14) 


=*./( w ) = 


A 3 8 k j + A 2 *{ 1 } + Attff + 


A 4 -A, 


(7.15) 


Here, S k ,i denotes the Kronecker delta; and the constants are given by 

*£} = r 4 n] (x)Mx)dx, n = 1,2,3, (7.16) . 

’ Jo 
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where here and in Eq. 7.14 t/»j[ n ^(x) = d n y?k{x)/dx n . 

Note that the terms described in Eqs. 7.14 and 7.15 consist of rational functions of 
the frequency u ; combined with pure delays. Thus, this stucture is straightforward to 
realize in hardware or software using linear, lumped-parameter networks with delays. 
Note that there are only two fixed delays, e ,u> ^ c and e~ lLj/Uc , independent of the choice 
of k and /, so that the entire stucture (At./fcu) + Aj k (u>) + E* i ;(u.’) + can be 

realized as four seperate (rational) matrix networks combined with scalar time delays. 
The remaining part of r(u>) is the term F(u), which is not rational. In particular, for 
the subsonic case, F is given by: 

FM = {p 1 )- t. (T.17) 

V n = l 

and for the supersonic case, F is given by: 

(-.is) 

"c .nasi 

Since this component of the spectrum is nonrational, it cannot be synthesized using 
a linear, lumped-parameter model. Fortunately, this part of the spectrum is not tied 
to the spatial structure; therefore, it can be synthesized as a scalar random process, 
which can then be fed into the spatial structure via the rational network derived from 
the remaining terms. The sums in Eqs. 7.17 and 7.18 can be realized by producing 
independent processes having the individual spectra of the summands, and then adding 
them together. Thus, the key to simulation of the nonrational part of this excitation is 
to synthesize stationary random sequences with the generic spectral shape 

<f>(us) = 7rae -0, H (7.19) 

To consider the generation of this process in discrete time, we first note that the 
autocorrelation function corresponding to Eq. 7.19 is given by 

^TTW -°° <r<co ' (720 > 

When sampled at time instants spaced A time units apart, the resulting discrete-time 
process {A*} will have autocorrelation sequence 

< 7 - 21 > 
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where oa = af A, which has the corresponding spectrum 


^2, e ,nw cosha&(w-t 

^ U - 12^ i + (kla±y ~ *** sink o A w 


(7.22) 

Such a process can be produced by filtering a white Gaussian sequence with a linear 
filter having gain \Jd2 (cu). If causality is not required, this can be accomplished simply 

by using the filter with transfer function \J 4>^(u). If, on the other hand, causality is 
desired, then it is necessary to factor d>$ into its causal and anti-causal spectral factors, 
<f>£ and 6 ^ , respectively, a task that must be performed numerically. 


8 Benchmark Examples 

In order to study the control problem in the above models, it is useful to specify bench- 
mark values for the various parameters arising in the preceding sections. 

For the subsonic case, we can choose among several scenarios studied in [19], We as- 
sume the same panel parameters as in the supersonic case, and fix the following param- 
eters for the boundary layer [19] (see Eq. 43): v = llOm/s, S' = 0.39cm, U c = 
0.8t>, and 6 = 17S'/U C . Again, we choose negligible damping (as in [19]) although 
this might perhaps be changed to a more realistic value. Within these benchmarks, we 
can consider the hierarchy of models listed in Table 8.1. The shaded region represents 
those models investigated in this Phase I report. 

In studying the beam, we will take intensive quantities to have the same numerical 
values as their higher-dimensional counterparts in the benchmark model. Note that 
the above hierarchy includes non-Gaussian excitation fields (such as found in a transi- 
tional boundary layer) as described in preceding sections. This hierarchy may also be 
extended to include nonlinear dynamics of the beam and panel by replacing the linear 
equations with appropriate nonlinear modifications (see, e.g., [7]). Consideration of the 
models in this order allows us to study the most basic phenomena first. In the simplest 
case, the primary phenomenon being studied is that of a second-order, linear stochastic 
differential equation, driven by a nonrational (scalar) random process with the model 
Eq. 43. By performing the analysis on a modal basis, we will see the effects of the 
temporal spectral shape on the control problem without complications due to spatial 
complexity. 

The latter will arise in the multimodal version of this same problem, and the com- 
plexity can be built up by adding modes until sufficiently rich modal structure is ob- ' 
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Table 8.1 Hierarchy of Models; Shaded Region Represents Those 
Models Investigated Experimentally in this Report 


Beam 

Subsonic ^ 

Gaussian . 1 111 

Single-Mode 
■ Multi -Mddi£; ' ; 

NOT-Gaus^art;l 

\ '' • ' 

Single-Mode: 

Supersonic 

Gaussian 

Single-Mode 

Multi-Mode 

Non-Gaussian 

Single-Mode 

Multi-Mode 

Plate 

Subsonic 

Gaussian 

Single-Mode 

Multi-Mode 

Non-Gaussian 

Single-Mode 

Multi-Mode 

Supersonic 

Gaussian 

Single-Mode 

Multi-Mode 

Non-Gaussian 

Single-Mode 

Multi-Mode 


tained. Consideration of the supersonic case within the beam model will require the 
development of a Green's function for the two-dimensional acoustic potential of Eqs. 
7.3-75 (analogous to that of Eq. 3.5). Work related to this problem can be found in [9]. 
Finding this Green's function is an intermediate task to finding a suitable 
representation for the acoustic damping in terms of the modal response of the beam. 
The study of the plate is a conceptually straightforward extension of the study of the 
beam. However, it is anticipated that the additional complexity will provide 
several challenges in implementing this extension. One important issue is the 
eigenfunction decomposition of the plate motion, which in the case of proportional 
damping must be accomplished by either numerical approximation, asymptotic 
expansion, or some other approximation method. One promising such method is that 
applied in [25], in which the modal behavior is approximated by separating the 
modal structure into one-dimensional (x and y) components. For non-proportional 
damping, eigenfunction decomposition is not possible. These are issues to be addressed 
in the follow-on to the current study. 

9 Experimental Results 

In this section we describe the results of computer simulation experiments that were 
conducted on a clamped beam as a first step toward adaptive nonlinear feedback 
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controls for reducing the vibratory response of aircraft panels subjected to a turbulent 
boundary layer. Investigation is limited to the shaded region in Table 8.1.* 


9.1 Control System Simulation 

We consider in this section the details of experiments at the first level in the 
hierarchy described in the preceding section. That is, we consider the control of the 
fundamental mode of a clamped beam subjected to a subsonic flow. 

A system model for this problem is the following forced, undamped scalar second- 
order ordinary linear differential equation (see Eq. 5.1): 


f PiYt 2 + VTt + DA >] W ^ = ^ + 


(9.1) 


where and D are the density and rigidity, respectively, of the beam; X, is the 
pfjncjpal eigenvalue of the biharmonic operator in one dimension (i.e., Eq. 7.6 with Jc 
= 1); and where VJ v f v and Cj denote the projections of the beam displacement, 
forcing field, and control applied, respectively, onto the principal spatial 
eigenfunction, \f/i, of the beam. Unless otherwise specified, Y ra = 0.5 to , where a> d = 


V * " is 1116 damped natural frequency of mode m, co m is the undamped natural 

frequency of mode m, and >s the damping ratio; these terms were discussed in 

Section 3, and are discussed further in Appendix A and B. 

Note that the homogeneous version of Eq. 9.1 is a damped oscillator with resonant 


frequency 


_ V PA./p! 

2n 


(9.2) 


For the benchmark 
eigenvalue is 


parameter values set forth in the above section, the principal 



(9.3) 


and thus the resonant frequency of Eq. 9.1 is 


Vj = 90.86 Hz. 


(9.4) 


The resonant frequencies are tabulated in Table 9.1 for the different modes. 
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Table 9.1: First Tsn Eigenvalues and Corresponding Natural Temporal 
Frequencies of Benchmark Beam 


Mode Number 

(m) 



G) m (rad/sec) 

1 

5.7842 

90.862 

570.90 

2 

43.948 

250.46 

1,573.6 

3 

168.92 

491.02 

3,085.1 

4 

461.58 

811.68 

5,099.9 

5 

1,030.0 

1,212.5 

7,618.3 

6 

2,009.3 1 

WESaam 

10,640 

7 

KE3ZSI 

2,254.7 

14,166 

8 

5,875.9 

2,896.0 

18,196 

9 

9,168.4 

3,617.5 

22,729 

10 

13,682 

4,419.1 

27,766 


To design an experiment for assessing the effectiveness of adaptive control 
according to Eq. 9.1, we first convert this model to discrete time. As will be seen, a 
sampling rate f s of 10,000 samples per second (i.e., a sampling interval of 100 psec.) 
will provide ample margin for implementing a digital control on Eq. 9.1; this is the 
sampling rate used in the simulations. On defining the sampled displacement and 
input sequences as 


o 

II 

■V 

< 

sT 

n 

3 

( 9 . 5 ) 

and 


f k =fi(k A), k = 0,1 , ... 

( 9 . 6 ) 

where A = l// s , and k is now the time index, we can write a discrete-time plant model 

corresponding to Eq. 9.1 as (see Appendix A for derivation) 


r Wjk+ll r * -aA [Y jin(o rf A)/pj - a sinidtjA) + ©</ cos({djA)]Io>j 


L w k + 1 J L - B}^e~ a ^sin(o>dA)/(o c i e 


[»£]* [a/“,] a *[4/°,] ci 

( 9 . 7 ) 
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for k = 0, 1, . . where itf* is the state variable representing the time derivative of 
Widt) at the sampling times t = kA. ... ' 

For the purposes of experimentation, we assume that we can observe directly the 
motion of the beam projected onto the fundamental modes (i.e., ideal sensing). The 
case of practical, discrete sensors is discussed further in Section 11.3. We also assume 
that we can measure the sequence W\, Wi, ■ . without error (i.e., the noiseless case). 
The situation in which there is measurement error is a straigh tforward modifica tion 
of this case. Moreover, we assume that we can apply a control C* at each sample time 
k. For a stochastic plant such as Eq. 9.7, the control sequence (c^j must, of course, be a 
function of the system outputs; i.e., for each k, C* will depend on previous W*s. 

Fig. 9.1 illustrates a block diagram of the complete simulated system with unity 
feedback control for each vibration mode; note that each mode to be controlled has a 
similar control loop. Simulation of the dynamics is direct from the plant and 
controller equations. The reference input is taken to be zero since this is the desired 
amount of modal deflection. The sequence [u> m (t)}, the modal displacement of the m th 
mode, is thus a direct measure of the error signal. The error signal drives the 
controller, which synthesizes a modal actuation, C m (t), that is summed with the 
turbulent modal force, f m (t) L that together drive the plant (i.e., act on the beam). The 
action taken by the controller should be such that the modal displacement, zy m (f), 
remains close to the reference input, which is taken to be zero. 



Figure 9.1: Block Diagram of Control System for Modal Response of 
Beam 



( 9 . 8 ) 


i m 

m 


w k> 


Jt=l 


Assuming as an objective the minimization of the average vibrational energy 
over m modes, an optimal control strategy can be specified for this problem; note that 
<■> denotes the expected value of the bracketed quantity. Generally, since F(co) is not 

rational but instead exponential, this control will be produced by an infinite impulse 
response (IIR) filtration of the measurement sequence. An HR filter structure has 
particular advantages for the problem of interest here, in that it is easily adapted; it 
is also easily generalized for the nonlinear case in which the linear predictor is 
replaced with a nonlinear polynomial neural network predictor. 

If the forcing input were a rational process, then the optimum filter would have a 
recursive implementation as a Kalman state predictor. However, for the case at hand, 
the system is not rational, and the optimal controller will not, in general, have such a 
simple implementation. As discussed below, however, since the forcing function is 
being modeled as an ARMA process, the optimum stochastic controller will take the 
form of an optimum state estimator (e.g., a Kalman filter). 

In its adaptive form, the objective of minimizing Eq. 9.8 is replaced with that of 
minimizing the raw squared average: 


m 


m k 
Jt=l 


( 9 . 9 ) 


The basic design issue is that of creating an adaptation algorithm to minimize the 
objective of Eq. 9.9 within the dynamical model (i.e., Eq. 9.7). The spatial issues have 
been removed by the projection of the spatial structure onto the principal eigenspace of 
the biharmonic operator. In this way, we reduce the problem to the above-described 
scalar problem in which we have the scalar measurement sequence {U> t ). Note that, in 
a real system, an approximation to such a measurement can be produced by combining 
the outputs of an array of sensors distributed on die beam. The outputs of these sensors 
can be combined at each sampling time in a way that mimics the projection of the 
displacement onto the eigenfuntion y/ m : 

a 

tvjt) = \\i/ m (x)w(x,t)dx. ( 9 . 10 ) 

o 
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Note that from the viewpoint of the controller, the plant is^unknown. This makes 
developing a controller difficult, since to determine its coefficients, one would ideally 
know the desired output of the controller in real time. The first approach that was 
attempted in designing a controller was to try to predict/^f) based on knowledge of 

both the plant and the previous measurements of {w k }, * nten fi° n °f setting 

C (f) = -/ (f), where C m (t) represents an estimate of C m (t); the sequence of previous 
measurements [W k ] must be used instead of \f k ] since the latter is not available. In 
designing such a controller one obtains for each mode an equation of the form 



it 

t\t-2 + "/f-llf-2^ 2 ' 

f=l 


(9.11) 


where [zo { } is the modal beam displacement, \f t \ is the modal forcing signal (i.e., 
boundary-layer pressure), and j M _ 2 and/ M , ,_ 2 are optimum two-step-ahead and one- 
step-ahead predictors of {/), respectively. 

When an FIR transversal filter based on this approach was implemented, it was 
found to be highly unstable; any error in predicting \f t \ introduced an error in the 

control signal input to the plant; this in turn increased the modal displacement error, 
which then had to be compensated for along with the turbulence force at the next time 

step. 

An alternative strategy was then implemented which proved to be highly 
successful. It is based on the notion of making no assumptions regarding the temporal 
relationship between c m (t) and The controller is implemented as an HR filter 

using a separate dynamic polynomial network nodal element, as shown in Fig. 9.2, to 
control the displacement in each mode. Note that an HR filter was used instead of an 
FIR filter since it was expected to reduce the number of terms necessary in the control 

system. , „ 

When designing an estimation neural network for the purpose of smoothing, 
filtering, or prediction, the structure of the network can be established using the 
Barron Associates, Inc. (BAI) Algorithm for Synthesis of Polynomial Networks 
(ASPN-11) [13]. This software package is used to structure and parameterize static 
(i.e., FIR) networks; similar algorithms exist for structuring and parameterizing 
dynamic (i.e., HR) polynomial networks [1]. However, because the modal 
displacement error input to the controller is a function of previous control outputs, such 
algorithms cannot be used to structure and parameterize the controller a priori-, 
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instead, the proper controller must be identified during the process of controlling. In 
practice, BAI has found that high-order controllers are generally not necessary. 
Accordingly, a simple dynamic nodal element feeding back the two previous control 

outputs, C m (f-1) and C m (t-2), and the three previous measurement values, 
and u; m (f-2), was attempted initially. Because this structure performed well, 
attempt was made to reduce the number of measurement value feedback terms. It was 
found that the minimum adequate controller is one that feeds back the two previous 
control outputs and the two previous sensor measurements. No attempt was made to 
reduce the number of previous control outputs that are fed back: two are the minimum 
necessary for constructing an oscillator which, in the case of no damping, captures the 
basic behavior of the plant. 



* shift register 

Figure 9.2: Dynamic Polynomial Neural Network Nodal Element; e(t) 

Represents the Modal Displacement Error and c(t) the 
Controller Output 

In all experiments, this simple controller was trained off-line using Gaussian 
white noise (GWN) as the input process f(t). GWN simulates any (broadband) 
continuous input disturbance to an arbitrary degree of accuracy. This implies that, for 
a plant with a given eigenstructure, in a mean-squared-error sense and without a priori 
knowledge of the temporal excitation function (other than its continuous as opposed to 
discrete nature), there is a single ideal control function. The spatial forcing function is 
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already taken into consideration in the modal response function as defined in Eq. 7.8. 
With this approach, the temporal excitation signal simulates any and all continuous 
flow conditions and thereby, indirectly, any and all boundary-layer excitations. 

With all parameters initially set to zero, the simulated beam was exposed to 0.25 
sec. of GWN, and the modal deflection was noted. The controller parameters were 
then adapted using a Guided Random Search (GRS) technique; the process is 
functionally similar to using the Least Mean Squares (LMS) procedure, but has been 
found to converge much more rapidly than the LMS algorithm. The simulated beam 
was then exposed to additional independent GWN sequences of the same duration, 
each time using a controller based on the previously adapted parameters to attenuate 
the modal response of the beam. Adaptation of the controller parameters continued in 
this fashion until the sum of the squared errors over each 0.25 sec. excitation interval 
had decreased to a steady-state level, which was generally between three and six 
orders of magnitude below the uncontrolled case (i.e., 27 to 60 dB vibration reduction). 

As discussed above, all mode controllers implemented in this study utilized only 
four parameters; two parameters were used to weight the previous control outputs, 
c m (t-l) andC m (t-2), and two to weight previous measurements of the modal 
displacement of the beam, w m (t-l) and w m (t-2). These controllers required 
approximately 1000 GWN sequences of length 0.25 sec. to achieve the above-stated 
error reduction, independently of the mode, m, being controlled. This suggests that 
real-time, off-line training of the network requires approximately four minutes. Note 
that training is done only once off-line to obtain nominal control parameters for each 
mode to be controlled; in addition, training of all mode controllers may be performed in 
parallel. Also, the coefficients determined in training the controller for each mode 
were used under all flow conditions. Because the controller is trained on GWN, it is 
general enough to handle essentially any broadband flow condition. 

Since the plant in the simulation is stationary, no on-line adaptation was 
required. An adaptive controller may be implemented, however, by continuously 
observing, on-line, the decrease or increase in the modal response of the beam to 
parameter changes, as was done during off-line training. In essence, the controller 
continually searches for new parameters to improve performance. In general, the 
search space is a nonlinear function of the parameters; thus, search techniques such as 
GRS may be used. Alternatively, Levenberg-Marquardt (LM) or LMS algorithms may 
also be used, based on the equations that describe the plant or on-line approximations, 
respectively. Since the plant equations are well known, or can be learned (e.g* 
estimated) inductively based on empirical data, an LM algorithm is a practical way 
to achieve significantly more rapid convergence of the parameters both on- and off- 
line. 
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9.2 Beam, Single Mode, Subsonic Flow 

9.2.1 Laminar Boundary Layer 

For laminar flow, the excitation process /(f) is assumed to be sinusoidal, as discussed in 
Section 4. Fig. 93a illustrates the response of the beam to pure laminar flow when no 
control is applied, where the sinusoidal excitation frequency was selected to 
correspond with the resonant frequency of the first mode (viz., 90.86 Hz). The 
amplitude of the excitation was scaled to equal that of the turbulent filter discussed in 
the next section. From the figure it can be seen that the projection of the beam 
displacement, XV v onto the first principal spatial eigenfunction, y/ l , grows in an 
oscillatory fashion until a steady state is achieved. Fig. 9.3b illustrates the modal 
response of the beam when control is applied continuously during beam excitation. 
Table 9.2 provides the average sum-square modal displacement error of the beam for 
the first mode, as computed over the last 0.1 sec. of data in the figures. From the table 
it can be seen that the controller reduces the average sum-squared error of the modal 
response of the beam by 60.3 dB. 



time (sec) 

Figure 9.3a: Uncontrolled First Mode Response of Beam to Laminar 
Excitation 
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Figure 9.3b: Controlled First Mode Response of Beam to Laminar 
Excitation; Note that Linewidth has been Increased to 
Make Response Visible 


Table 9.2: Average Sum-Squared Modal Displacement for Mode One 

with and without Control^ 


Type of 
Excitation 

Average Sum- 
Squared Error 
(without 
Control) 

Average Sum- 
Squared Error 
(with Control) 

Vibration 
Reduction with 
Control (dB) 

Laminar 

2.43134 x 10' 10 

2.24831 x 10‘ 16 

60.3 

Turbulent 

7.25780 x 10 14 

6.70857 x 10' 17 

30.3 

Transitional ! 

4.87394 x 10' 11 

1.64366 x 10' 16 

54.7 


9.2.2 Turbulent Boundary Layer 

In exciting the beam with turbulent Bow, a major challenge is to simulate the forcing 
process. This part of die simulation requires generation of random sequences with the 
power spectral density of the sampled continuous-time disturbance. For die first mode 
this spectral density is given by 

<Pd(co) = f s I - \ A (f s 0)), -1t<(D<1t, (9.12) 
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where Tj j is given in Eq. 7.13. Such a sequence can be generated by filtering a white 
Gaussian process with a filter whose transfer function is 


H(co) = yj (Pdico), -K<ca<n. (9.13) 

This filtering can be performed in the frequency domain, as shown in Fig. 9.4, by first 
generating a white Gaussian sequence; applying a fast-Fourier transform (FFT); 
multiplying the transformed sequence by H(<y); and then retransforming back to the 
time domain. Implementation of this model directly in the time domain is also 
straightforward, as will be seen. 



Figure 9.4: Simulation of Turbulent Forcing Function 


From Eqs. 7.10 - 7.18, the turbulent forcing function is seen to be a Gaussian 
random process with power spectral density: 


ri,i(<a) = 2 F(co) [Re{Ai t i} + Refs^}], (9.14) 


where Aij, S\ t \, and F, are as defined in Eqs. 7.14, 7.15, and 7.17, respectively. The 
function F(dtf) is plotted versus CO in Fig. 9.5 for the subsonic benchmark parameters. 
Note that this part of the spectrum is relatively broadband, having a (one-sided) 
half-power bandwidth of approximately 3000 rad./sec. (500 Hz). For reference, the 
fundamental frequency of the unforced oscillator occurs at 570.9 rad./sec., as shown in 
Table 9.1. 

The remaining term in Eq. 9.14 is a rational function of frequency. Its primary 
behavior as a function of frequency can be predicted from the function 

1 

I A 4 - Ai I 

as seen in Eqs. 7.14 and 7.15. 


(9.15) 
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Figure 9 . 5 ; F((o), the “Irrational" Part of the Forcing Spectrum for 
Subsonic Benchmark Parameters: 5* = 0.0039m, v = 

11 Om/s, A', = 0.24, A’ 2 = 1.08, A 2 = 1.80, K\ = 0.47, K 2 = 3.0, 

and JCj- 14.0 

Fig. 9.6 shows the behavior of the real part of the rational term Ai.i for the 
subsonic benchmark parameters... Note that this function has a much narrower 
bandwidth than that of F; thus, it will determine the overall bandwidth of the 
forcing function. Fig. 9.7 illustrates the behavior of the real part of the rational term 



Figure 9.6: Real Part of A t ^0)), a Component of the “Rational" Part of 
Forcing Spectrum for Subsonic Benchmark Parameters: 

A = g + iw/U c ,g = l/(ir c 0), li^e = 178*, u e = 88 m/s, and 

5* = 0.0039m. 
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0.1 



© (rad /sec) 

Figure 9.7: Real Part of £ 1,1 (©), a Component of the “Rational” Part of 
Forcing Spectrum for Subsonic Benchmark Parameters: 

A = g + iiv/U c ,g = 1 J(U C Q), U C Q = 175*, U c = 88 m/s, and 
5* = 0.0039m. 


<21,1/ again using the subsonic benchmark parameters. The complete turbulent 
excitation forcing function of Eq. 9.14, which is real everywhere for mode one, is 
plotted in Fig. 9.8. 



<d (rad/sec) 

Figure 9.8: Spectral Density of Turbulent Excitation for Mode One 
Based on Subsonic Benchmark Parameters 
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It should be noted that, in view of the fact that the function F(co) is relatively 
broadband, a reasonable approximation in computing T would be to assume that 
F(co) is, in fact, constant. This approximation models the forcing process as a rational 
system driven by white noise. In particular, under this assumption, the spectrum 
F i,i(<y) becomes a ratio of polynomials in Cifi. These polynomials can be factored into 
causal and anti<ausal parts, the causal part of which defines a filter generating the 
desired forcing process when driven by white Gaussian noise, (see, e.g.. Chapter 5 of 
[16].) This was not the approach taken herein, however. As discussed earlier in this 
section, Eqs. 7.14, 7.15, and 7.17 were solved analytically for the first two modes. 
These solutions are rather extensive and are therefore not given here; they were 
generated using Mathematica, a software package for performing symbolic (as 
opposed to numeric) calculations by computer. Instead, their spectra are given 
graphically: for the first mode, there is a single filter term Ti^ca), which was 
depicted in Fig. 9.8. This filter was numerically simula ted f or — 7C < CO < n, and the 
inverse Fourier transform computed using an FFT algorithm, to obtain the turbulent 
excitation impulse response function in the time domain. A white Gaussian sequence 
was then filtered in the time domain using an FIR filter with this impulse response, to 
produce the turbulent excitation signal. Fig. 9.9 shows a typical result in exciting the 
panel with turbulent flow. Note that in the uncontrolled case (Fig. 9.9a), the modal 
response of the beam achieves a steady state having significantly less modal 
displacement than was true for the case of laminar flow. When control is continuously 
applied (Fig. 9.9b), the average sum-squared modal displacement is reduced by 30.3 
dB. Fig. 9. 10 illustrates the power spectra corresponding to the first mode 
displacement depicted in Fig. 9.9. Note that the controller attenuates significantly 
vibration at all frequencies, not only the resonant frequencies, of the beam. 

9.2.3 Laminar-to-Turbulent Boundary Layer Transitions 

The simulation of the excitation function in transitional flow was performed using a 
binary gating field that switched intermittently between the two equilibrium 
boundary layers (i.e., pure laminar flow and pure turbulent flow). The binary gating 
field was modeled as Markovian in time, with transition statistics depending on the 
probability of being in either state. For these experiments the probability parameter, 
p, was set equal to 0.5, indicating equal likelihood of laminar and turbulent flow 
conditions. Setting p = 0.5 produces the most intermittent flow regime possible with 
such a simple model. 

The result of intermittent flow gating in the uncontrolled case is shown in Fig. 
9.11a. When control is applied the response is as given in Fig. 9.11b, where it is seen 
that the resultant first mode response is reduced by 54.7 dB. 
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Figure 9.9a: Power Spectrum of Modal Displacement with Turbulent 
Boundary Layer for Mode One without Control 
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Figure 9.9b: Power Spectrum of Modal Displacement with Turbulent 
Boundary Layer for Mode One with Control 
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Figure 9.10a: Uncontrolled First Mode Response of Beam to Turbulent 
Excitation 



Figure 9.10b: Controlled First Mode Response of Beam to Turbulent 
Excitation 



Figure 9.11a: Uncontrolled First Mode Response of Beam to Laminar 
to-Turbulent Excitation 



Figure 9.11b: Controlled First Mode Response of Beam to Laminar-to- 
Turbulent Excitation 
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As discussed earlier, laminar-to-turbulent boundary-layer transitions represent a 
non-Gaussian excitation which, in general, will require the use of nonlinear elements 
in the adaptive PNN controller. In our simulations, however, it was found that linear 
networks adequately attenuated the modal displacement We believe that this was 
the case because each mode of the beam represents a narrowband filter with respect to 
the broadband turbulent excitation. The bandwidth of a system is inversely 
proportional to its time constant, and this plant has a relatively Targe time constant* 
The consequence is Sat the effects of the input excitation are averaged in time and 
therefore tend to a Gaussian distribution, for which case a linear controller is optimal. 

9.3 Beam, Multi-Mode, Subsonic Flow 

The introduction of multiple modes complicates the above problem in that the scalar 
Eq- 9.1 is replaced by the vector Eq. 5.1. The turbulent forcing field then is a vector 
Gaussian process with a spectral density matrix described as in Eqs. 7.14, 7.15, and 
7.17. The control is now a vector process and the adaptation becomes a vector 
adaptation. For the various flow conditions illustrated below, results are provided for 
the second mode only since the vibration attenuation in each mode is essentially 
independent of the other modes; therefore, results obtained for the second mode were 
similar to those demonstrated above for the case of single-mode control. 

9.3.1 Laminar Boundary Layer 

Laminar flow was again simulated as a sinusoidal process, this time with a frequency 
equal to that of the resonant frequency of the second mode of the beam (viz., 250.46 
Hz). Fig. 9.12a illustrates the response of the second mode of the beam to laminar flow 
when no control is applied. The projection of the beam displacement, W 2 , onto the 

second principal spatial eigenfunction, yf 2 , of the beam grows until it achieves a 
steady state! The response of the beam when control is applied is shown in Fig. 9.12b. 
From Table 9.3 it can be seen that the average vibration reduction as computed over 
the last 0.1 sec. of the data was reduced by 51.0 dB. 


9.3.2 Turbulent Boundary Layer 


To simulate turbulent flow, Eqs. 7.14, 7.15, and 7.17 were solved analytically for T. For 
the second mode, T is a matrix of dimension 2x2 that cross-couples the random 
processes exciting each mode independently; specifically, the boundary-layer forces 
are generated using the following equation: 



fuirCWN! 

1 2,2 JLgwn 2 


} 


(9.16) 
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Figure 9.12a: Uncontrolled Second Mode Response of Beam to 
Laminar Excitation 



Figure 9.12b: Controlled Second Mode Response of Beam to Laminar 
Excitation; Note that Linewidth has been Increased to 
Make Response Visible 
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Table 9.3: Average Sum-Sq uared Moda l Displacement for Mode Two 

with and without Control ; . 


Type of 
Excitation 

Hi 

Average Sum- 
Squared Error 
(with Control) 

Vibration 
Reduction with 
Control (dB) 


■ESSEH 

1.29531 x 10' 16 

51.0 

Turbulent 

8.25576 x 10‘ 14 

6.62533 x 10' 17 

31.0 

Transitional 

■mss 

1.34811 x 10' 16 

45.4 


where GWN l and G WN 2 are independent Gaussian white noise processes. Note that 
in the unimodal case, GWN 2 = 0, so that only the F t j term was needed. In the 
multimodal case, creating the turbulent excitation amounts to filtering independent 
GWN sequences using the multi-dimensional filter T. 

For the second mode there are four terms in die spectral density matrix; in general, 
for m modes, the matrix has m 2 terms, growing as die square of the number of modes. 
As m increases, the analysis and computational effort to generate die turbulent forcing 
functions having the appropriate spectral characteristics grows dramatically. For 
this reason, the number of modes studied in this Phase I report was limited to die first 
two. 

As mentioned earlier, the individual terms of T, namely r j j, T l2 , r 21 , and 
r, ,, were generated based on Eqs. 7.14, 7.15, and 7.17. The spectral density F j j was 
discussed earlier in Section 9.2.2 in describing the turbulent excitation for mode one; 
recall that Tj { in Eq. 9.14 is real because the imaginary terms cancel in Eq. 7.13; this 
is guaranteed in this case because A— = Aj j and Sj j = «3jj. For the same reason, T 2 2 
will also be real; Fig. 9.13 shows a graph of T 2 2 (£0) vs. <0. The remaining terms in 
the turbulent filter matrix T are F, 2 and T 2 j, which are complex since the 
imaginary terms do not cancel, because in general A- * A- and is- * j indeed, 
Fj 2 and r 2 j are complex conjugates of one another, so that their imaginary parts add 

rather than cancel. Graphs of the real and imaginary parts of T 12 and F 21 are given 
in Figs. 9.14 and 9.15 respectively. N ote from Eq. 9.16 that the forcing vector / is 
complex. Because we require a real excitation sequence, only die real parts of T 12 and 
T- were used to generate the off-diagonal elements in the turbulent excitation 
matrix. The T matrix was simulated as described earlier to create time-domain filters 
that were applied to the Gaussian white noise sequences. 
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Fig. 9.16a illustrates the case of no control for the second mode where the beam 
was excited using the simulated turbulent flow process defined by T. The 
corresponding controlled case is shown in Fig. 9.16b, where the vibration is seen to 
have been attenuated by 31.0 dB. The corresponding natural and controlled vibration 
spectra for this mode are shown in Fig. 9.17. 



to (rad/sec) 

Figure 9.13: Spectral Density Matrix Term T 2 2 








time (sec) 

Figure 9.16a: Uncontrolled Second Mode Response of Beam to 

Turbulent Excitation 



Figure 9.16b: Controlled Second Mode Response of Beam to Turbulent 

Excitation 
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Figure 9.17a: Power Spectrum of Modal Displacement with Turbulent 
Boundary Layer for Mode Two without Control 
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Figure 9.17b: Power Spectrum of Modal Displacement with Turbulent 
Boundary Layer for Mode Two with Control 
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9.3.3 Laminar-to-Turbulent Boundary Layer Transitions 
To simulate laminar-to-turbulent boundary layer transitions, the same binary gating 
field discussed in Section 9.2.3 was used. The response of the second mode of the beam 
with no control is illustrated in Fig. 9.18a. The controlled case is shown in Fig. 9.18b, 
where it can be seen that the vibration response has been reduced by 45.4 dB. 

9.4 Vibration Suppression as a Function of Damping 

To study the reduction in modal displacement in turbulent flow as a function of 
damping, the damping ratio, C used in the plant model was varied from 0.001 to 0.1. 
Typically, the range of damping ratios that one would expect to see in panels is on the 
order of 0.001 to 0.01. The controller polynomial neural network was trained using the 
power-law proportional damping relationship that is discussed in Section 3 and in 
Appendices A and B; the trained controller was then tested on plants having different 
levels of damping. As demonstrated in Appendix B, power-law proportionate 
damping provides an effective damping ratio of = 0.007265 at the first-mode 
resonant frequency, and - 0.003696 at the second-mode resonant frequency. Figs. 
9.19a and 9.19b illustrate, for the first two modes respectively, that the reduction in 
vibrational energy due to active control decreases with increasing damping ratio, as 
passive damping alone reduces the amplitude of vibration to be controlled. 

It is interesting to note that even in the case of nolJamping, the controller is still 
able to keep vibration to a negligible level, si mila r to t ha t demonstrated in Sections 
9.2 and 9.3. In the undamped case without control, the modal displacement continues 
to grow without bound for the case of intermittent flow with p = 0.5, as shown in Figs. 
20a and 21a for the first and second modes respectively; with control, the modal 
displacement remains negligible, as illustrated in Figs. 20b and 21b. As a result, the 
relative reduction in vibration level increases without bound as damping goes to zero. 
Note that the controller in the above undamped case was trained using power-law 
proportionate damping; it was not retrained for the undamped plant. Similar results 
were obtained for the cases of pure laminar and pure turbulent flow. 
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Figure 9.18a: Uncontrolled Second Mode Response of Beam to 

Laminar-to-Turbulent Transitional Flow 



Figure 9.18b: Controlled Second Mode Response of Beam to Laminar- 

to-Turbulent Transitional Flow 
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Figure 9.19a: Reduction in Average Vibrational Energy in First Mode 

with Turbulent Flow 



Figure 9.19b: Reduction in Average Vibrational Energy in Second Mode 

with Turbulent Flow 
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time (sec) 

Figure 9.20a. Uncontrolled First Mode Response of Beam to Laminar 
to-Turbulent Transitional Excitation with Zero Damping 



time (sec) 

Figure 9.20b. Controlled First Mode Response of Beam to Laminar-to' 
Turbulent Transitional Excitation with Zero Damping; 
Note that Linewidth has been Increased to Make 
Response Visible 
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time (sec) 

Figure 9.21a: Uncontrolled Second Mode Response of Beam to 

Laminar-to-Turbulent Transitional Excitation with Zero 
Damping 



Figure 9.21b: Controlled Second Mode Response of Beam to Laminar- 

to-Turbulent Transitional Excitation with Zero Damping; 
Note that Linewidth has been Increased to Make 
Response Visible 
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10 Conclusions 

Adaptive nonlinear polynomial neural networks were used to control vibration of a 
simulated aircraft structure (viz., beam) resulting from acoustic excitation by the fluid 
boundary layer. The objective was to control actively the coupling of sound from the 
fluid to the structure. The acoustic excitation in this project included laminar, 
turbulent, and laminar-to-turbulent transitional flows, representing respectively, 
narrowband-Gaussian, broadband-Gaussian, and broadband-non-Gaussian excitation. 
Excellent attenuation of 30 to 60 dB was achieved, even under conditions of no natural 
damping. 

11 Future Extensions 

11.1 Theoretical Analysis 

Due to the necessarily limited scope of the Phase I study, the concept demonstration 
was restricted to a simpler, lower-dimensional framework than can be reliably used to 
model actual aircraft panels. To bring the promise of active vibration control into 
practice, it is necessary to continue this study using more realistic panel models, and to 
verify the predicted control capabilities through experiment. The following are 
among the topics to be investigated further: 

11.1.1 Plate (versus Beam) Models 

The Phase I effort focused on the problem of a vibrating beam rather than a vibrating 
plate, the latter of which more accurately models aircraft panels. Thus, to bring these 
techniques to the point of applicability, it is necessary to add an additional 
(transverse) dimension to the structure being controlled. There are three primary 
issues that will arise because of this increase in dimensionality. One issue is that, 
unlike the case of a clamped beam, the eigenstructure of a clamped plate is not known 
in terms of readily calculated functions. This is primarily an analytical problem that 
can be overcome using approximation techniques. Another, more substantive, issue is 
that there may be effects in the control of a vibrating plate that do not appear in a 
vibrating beam. A further issue is that the possible configurations of discrete 
actuators and sensors (discussed below) will be quite large; thus, considerable 
investigation will be necessary to select optimal geometries. 
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11.1.2 Panel Nonlinearity 

The linear plate equation only approximates the true dynamics of a vibrating plate. 
This approximation is accurate for small deflections and linear-elastic materials 
where shear distortions and rotary inertia are negligibly small. Thus, in the problem 
of interest here, this approximation should be accurate for vibrations that cause small 
panel displacements. If the plate vibration is regulated effectively, a linear model 
will provide an accurate description of nominal plate behavior. On the other hand, 
even in a regulated plate, there is the potential for larger plate excursions if the 
excitation field is sufficiently intense. Such excursions will drive the plate into 
nonlinear behavior, and thus it is of interest to model this nonlinearity and to consider 
its effects on controller design and performance. The mathematical modeling of 
nonlinear plate behavior has been described in several works, including [3, 7]. There 
are several potential sources of nonlinearity in the plate behavior. Perhaps the most 
important of these phenomena for the thin, clamped panels often used in aircraft is in- 
plane stretching, which causes additional stiffness when the panel deflection is 
significant relative to the plate thickness. This stiffness increases with 
displacement, thereby introducing nonlinearity. For sufficiently large displacements 
(relative to plate thickness), this membrane action can actually dominate the bending 
action described by the linear plate equation. This phenomenon can be modeled 
mathematically using classical equations developed by von Karman [3]. In particular, 
the plate equation for larger deflections is obtained by adding a bilinear term to the 
right-hand-side of the linear plate Eq. 2.1. This term is given by: 


5 P 


"3 2 F 3 2 w 3 2 F d 2 w 

3 y 2 3x 2 + 3x 2 3 y 2 


3 2 F d 2 w' 
^ 3x3 y 3x3 y ' 


( 11 . 1 ) 


where 8p is the thickness of the plate, and F is Airy's stress function, which is 
quadratically related to the deflection w through the PDE: 



d 2 xv 3 2 zv 
'dx 2 3 y 2 } 


( 11 . 2 ) 


with E denoting Young's modulus. 

Subsequent efforts should consider these and related nonlinear phenomena both in 
terms of their effects on the linear control of vibrating panels, and in terms of their use 
in the design of appropriate nonlinear active vibration control laws for situations in 
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which the linear designs fail to exhibit robustness with respect to this nonlinear 
behavior. For weak nonlinearity, perturbations on the modal approach are suitable 
for analysis of these problems. However, to simulate accurately the behavior of a 
strongly nonlinear plate, it may be necessary to abandon the eigenfunction methods 
that are useful in the linear model and go instead to finite-difference methods of 
propagating the corresponding nonlinear PDEs. Such simulation methods will also be 
useful in other aspects of the follow-on studies. 


11.1.3 Accurate Sensor/ Actuator Models 

Practical sensors (e.g., polyvinylidene fluoride (PVDF) films) measure strain rather 
than displacement, and the control effects produced by practical actuators (e.g., lead- 
zirconate-titanate (PZT) ceramics) are bending moments rather than direct 
counterforces to the boundary layer. Although it is possible to produce sensors and 
actuators that couple into the fixed spatial modes of a plate (see, for example, [2, 6, 
10, 14, 23]), a more likely situation is that discrete sensors and actuators will be used. 
To increase the realism of analytical phases of control design and performance 
prediction, it is necessary to consider the constraints imposed by such physical sensors 
and actuators. Use of multiple sensors and actuators represents a first step toward 
systems that can be practically implemented. 

As a more realistic measurement model, we consider here the sensing/actuation of 
the field by discrete sensors /actuators having non-negligible dimensions, placed on 
the surface of the beam. (For the purposes of discussion, we consider a beam model.) 
One immediate effect of such a model is that, unlike the distributed case, it will not be 
possible to isolate modes in this context, and thus modal coupling through the 
controller must be considered. To model this situation, it is convenient to represent the 
dynamics of the beam in terms of an infinite-dimensional state-space model: 

4 - AX, + H (113) 

where the state vector X, is given by 


* 


w x {t) 

w 2 {t) 

w 2 (t) 
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with W k (t) and W k {t) denoting the k? h mode of the displacement field and its time 

a 

derivative, respectively, i.e., = J xv(x,t) Xf/^x) dx. A denotes the block- 


J 


diagonal matrix with 2 x 2 diagonal block 

r ° 1 

a -[- 

with Q) k = \ D \/Py As before, D denotes the beam stiffness, p ( its mass density 

a 4 

per unit length, \ k is the k 1 ^ eigenvalue of the operator -^ 4 , ^ k the damping ratio for 
the k! h mode, and/, denotes the vector whose k th component is the projection of the 
boundary layer pressure field onto the k th mode of the beam (i.e.,/ A (f) = Jg f(x,i) 
yr k (x) dx); B is the matrix 

A 
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Note that the modes are a fortiori uncoupled in tms moaei; ana, wun me aoove 
model based on perfect sensing and actuation, the modes remain uncoupled. On the 
other hand, if the beam motion is observed and actuated through practical sensors and 
actuators, mode-coupling will result. 

Suppose the beam motion is observed through a sensor occupying die region {x Q ^ X 

£ Xj} with 0 < x Q < *! < fl. The resulting sensor electrical signal r(t) equals K 

s(x,t)dx, where s(x,t) is the strain in the beam and K is a constant incorporating the 
efficiency of the sensor in converting strain into electrical energy. The strain is related 
through the displacement by 


, d 2 w(x,t) 
s(x,t) = K —^r 1 


(11.4) 
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where K' is a constant determined by the properties of the panel material. Thus, the 
sensor signal is given by 

r(f) = K s [w x (x v t) - U) x (x 0 ,t)], (11.5) 


where W x = 


d w(x,t) 

dx 


and K g = K K . 


Since iv(x,t) = Xjfc=i W k (t)\f/ k (x), the sensor signal can be written in terms 
of the state X f as 


r(t) = h. X t 


( 11 . 6 ) 


where h = K $ [(Y^ (*i) “ Vljc ( x o))' ^ ( x l) “ Vl,x •••! 

with y/fe - y k . Obviously, if we have multiple sensors, each one will satisfy 
such an equation. Thus, with L sensors we have a measurement equation: 


r 2 (t) 

<rtW> 


HX ( 


(11-7) 


with 

ft) 


\hJ 


( 11 . 8 ) 


and tty denoting the I th version of h. above. 

Similarly, practical actuation can be modeled as follows. An electrical signal into 
an actuator occupying foe region i/q < x ^ 1/j will produce a bending moment whose 
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magnitude is directly proportional to the voltage applied to the actuator; i.e., M(f) 
_ K a c(t), where c(t) is the control voltage across the actuator leads. (Formulas 
relating the constant K a to the material properties and geometry of the actuator are 

available.) This moment can be modeled as consisting of two equal-magnitude bending 
moments around the endpoints, y 0 and y v of the actuator. These bending moments can 

be taken to act in opposite directions so as to bend the center of the actuator down. In 
terms of forces driving the beam equation, these moments can be represented as four 

M(t) ' . 

point forces, each having magnitude ’ located at the points X - I/q - Ax, X - y 0 

+ Ax, x = y-i - Ax, and x = y 1 + Ax, where Ax > 0 is small relative to the sensor 
dimension. The forces at X = y 0 - Ax and X = y l + Ax are upward, and the forces at 
X = y 0 + Ax and x = - Ax are downward. Thus, the force created by the control 

signal c(f) applied to the actuator becomes 

fJx,t) = : ^Mi-y 0 -Ax)-^-y 0 + to) + 5(x-y 1 + Ax)- 
Ja 2Ax 

5(x-yi-Ax)], (119) 

where <5( • ) denotes the Dirac delta function. 

As an input to the beam state equation, this force translates into a vector j^(f) 
with k th component (i.e., input to the mode) 

r K fl c(f) 

/,,»(') - J /« WW *) dx - ~^r + **) - VkO/o - **) + 

0 

n(yi~ Ax) - + Ax)]. 

(11.10) 

The moment arm Ax can be eliminated from the above by considering the limit as 
it vanishes. In this limit (Ax -> 0), we have 

/„,*(<) = MO l^(y 0 ) - < nll > 
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Note that these forces enter the state equation as a term additive with the 
boundary layer vector With multiple non-overlapping actuators, we have the 

superposition of the signals from the various actuators entering the plant; i.e., 

M 

L«>= S/V«), (11.12) 

m=l 


where M is the number of actuators, and^^(f) is the actuation vector (as above) for 
the m ,h actuator. 

The above equations give us a complete parametric model for the transfer between 
the beam motion and die sensor outputs, and between die control signals and the beam 
motion. The control algorithm is mserted between these sensor output signals and 
control input signals. Thus, from these equations, the control synthesis and simulation 
can proceed. 

11.1.4 Supersonic Boundary Layers 

A significant increase in complexity results from the assumption of a supersonic 
boundary layer. Here, the damping effects of the acoustic potential must be calculated 
in modal form, a task that requires a calculation on die order of Eq. 7.10 but with 
singular integrands. 

11.1.5 Transitional Boundary Layers 

Transitional flow can be modeled as a binary gating field that switches intermittendy 
between die regimes of pure laminar and pure turbulent flow. In Phase I the binary 
field was modeled as being Markovian in time; follow-on investigations will require 
that this be extended to be Markovian in both time and space, which can be modeled 
as a dynamic Gibbs field. Transitional flow is expected to be a significant phenomenon 
only in subsonic flow regimes; in supersonic flow, the boundary layer is sufficientiy 
well established so as not to be intermittent. 

11.1.6 Multiple-Input, Multiple-Output Nonlinear Feedback Control 

Digital signal processors (DSPs) can be used to control the panel based on adaptive 
polynomial neural network algorithms. In the Phase I effort, the modal control 
algorithms used might be characterized as "selfish," in that the controller for each 
mode was concerned only with minimizing the modal displacement of that mode. In 
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pract ice, howeve r, effecting contr ol of one mo de will generally excite other modes, 
thereby causing spillover and increasing the displacement in these modes. It is 
desirable to have a more "selfless" control approach, one in which actuators work 
together as well as individually for the combined minimization of beam or panel 
modal displacement. 

Although individual control loops could be constructed for each mode over which 
control is to be exercised, a single multivariate controller is more powerful than 
multiple univariate controllers, as cross-couplings T>etween modes can be taken into 
consideration, as shown in Fig. 11.1; this represents a perfectly general multi-input, 
multi-output controller (MIMO) for the panel vibration. Note that there is no 
requirement that the number of force actuators be equal to the number of displacement 
sensors. - 



Figure 11.1: General Multi-Input, Multi-Output, Polynomial Neural 
Network (PNN) Mapping a V ector of L Sensor 
Measurements, x(t), Into a Vector of M Controllers, y.(t); 
Note that the Dynamic PNNs (see Fig. 9.2) may Contain 
Output Feedbacks and Internal Time Delays 


54 




Prior to the advent of neural networks, no general methods were available for the 
synthesis of MIMO, nonlinear, adaptive controllers. Neural networks may now be used 
to model implicitly the complex dynamic relationships and couplings between modes 
for completely general MIMO control. Recently [21], Barron Associates, Inc. (BAI) 
demonstrated that polynomial neural networks subsume Wiener-Volterra system 
identification techniques, which are often used for the identification of MIMO 
systems [20]; with PNNs, nonlinear and nonstationary (time-varying) systems, even 
those having infinite memory (e.g., oscillatory systems), may be readily identified. 
MIMO control represents an important research area having many potential 
applications. 

BAI believes that successful development of such algorithms is feasible today. 
For the control of panel vibration, a MIMO polynomial neural network will be 
realized by adapting the network to achieve a globally optimum solution. This 
solution approach was not attempted in the Phase I effort due to time and budget 
limitations. 

11.2 Laboratory Experiments 

Laboratory experiments might be designed to simulate the turbulent boundary layer 
excitation and control of a section of an airplane panel consisting of a rectangular 
frame covered with a thin aluminum alloy skin. Three sets of experiments are 
suggested, beginning with applications of the control method to simple beam structures 
subject to controlled random loadings (e.g., random, progressive plane waves) and 
progressing to high velocity turbulent boundary layer excitation of panels. 

11.2.1 Proof of Concept Experiments 

The initial series of experiments are designed to demonstrate that a polynomial 
neural network can control the dynamic response of a simple, multimodal structure 
subject to random loadings. The structure will consist of a beam (clamped /damped 
boundary conditions) mounted in one wall of an acoustic wave guide. The acoustic 
sources for the wave guide will consist of a series of horn drivers that generate plane, 
progressive waves with controlled amplitude and frequency spectra. The wave guide 
will be terminated anechoically to minimize end reflections. A series of sensors and 
actuators will be mounted along the outside surface of the beam in sufficient numbers to 
provide sensing and control in bandwidths containing high modal densities of the 
beam. It is suggested that approximately three sensor /actuator pairs should be used in 
this study. 
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Amplifiers for the actuators and signal conditioners for the sensors should be 
selected to be compatible with analog-to-digital (A/D) and digital-to-analog (D/A) 
computer control. The controller unit might consist of a digital signal processing (DSP) 
board interfaced with the A/D and D/ A converters. For given excitation conditions, 
die controller unit drives the actuators independently, with a phase and magnitude 
prescribed by the neural network algorithm. Parameters to be varied in the 
demonstration experiment include the amplitude, bandwidth, and bandwidth center 
frequency of the progressive acoustic wave field. For this particular excitation, the 
controller unit will drive the actuators to minimize the mean square velocity of the 
beam response at a sufficient number of sensor locations to produce a global minimum. 

11.2.2 Subsonic Turbulent Boundary Layer Excitation of a Uniform 
Beam 

In this series of experiments, the acoustic wave guide is fitted with a centrifugal 
suction fan to provide turbulent boundary-layer excitation of beam modes up to flow 
velocities corresponding to Mach 0.1. Hot-wire measurements are taken in the 
turbulent boundary layer to determine die spectral content and spatial structure of the 
flow. A measure of success of the adaptive controller will be its ability to reduce the 
response of the beam globally as a function of flow velocity and turbulence level. As a 
further test of the robustness of the controller, turbulent bursts will be introduced 
artificially upstream of the beam to provide potentially destabilizing fluid dynamic 
forces. 

For the experiments outlined in Sections 11.2.1 and 11.2.2, the controller strategy 
should be examined to determine its ability to control the response of a multimodal 
structure with high modal densities and subject to distri buted loading based on 
stationary, as well as time-dependent, spectra. Additionally, as a major part of the 
above effort, the success of the control strategy as a function of the number of sensors 
and actuators should be investigated. Particular attention should be given to setting 
guidelines for the number of actuators and sensors required in the frequency band 
containing a specified number of structural modes. 


11.2.3 High Velocity, Turbulent Boundary Layer Excitation of Panels 

Following the successful completion of the above experiments, a third set of 
experiments should be conducted to control the vibration response of a rectangular 
plate subject to simulated, turbulent boundary layer excitations at high flow 
velocities. This excitation should simulate the distributed and random loading 
conditions of a high velocity, turbulent boundary layer. Based on the results of the 
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experiments outlined in Sections 11.2.1 and 11.2.2, the control algorithms should then 
be extended to accommodate a larger number of sensor /actuator elements. For the plate 
experiments, it is suggested that approximately 100 sensor /actuator pairs be involved, 
forming a ten by ten grid. 
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Appendix A Derivation of Plant State Equations 

The vector second-order ordinary differential equation for the damped beam is: 

POP + TJi + DA ] ^ = M + ( A -l) 

where D is the plate rigidity; p j is the mass density (per unit area) of the panel; yis 
the damping law; A is the the diagonal matrix with diagonal elements representing 
the eigenvalues of the solution to V 4 ty/ = X\f/; W is the modal displacement; and /and 
C are the vector boundary layer excitation and control forces, respectively. 

Letting 

Xj(f) = w,(f), 

x 2 (t) = w l (t)=x l (t), 

* 3 (0 = 

X 4 (t) = W 2 (t) =X 3 (t), 
etc., 

the state equations for this system can be written in the form 

2£ = A £+B/+E£ ( A.2) 
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in which the A matrix is seen to be block diagonal with diagonal terms for the m th 
mode of the form 


r ° 

1 1 

‘ 0 1 

L ~ DX m'P\ 

-Y/Pi J ' 

-2C (0 m 


To put this model into the form 


z* + 1= hx* + /i + £1 


(A.4) 


for simulation requires the state transition matrix (j^, which is defined by 


<p k = £ _1 [(sl _ A)’ 1 ] t 


(A.5) 


where L is 
interval, I is 
gives 


the Laplace operator, s is the complex frequency, A is the sampling 
the identity matrix, and A is defined as above. For the m ,h mode Eq. A.5 

-1 


K - rl 


s + y/pj 


s 2 + ys /pj+ DX m /p l s 2 + ys/p 1+ DX m /p l 

zEhn/. El - 5 

L s 2 + 7 S/P 1 + DX Jp x s 2 + Tfs/ Pi+ DX m / Pl J 


(A.6) 


t =A 


Comparing the polynomial in the denominator of each of the above terms with the 
characteristic second-order equation S 2 + 2£ + CD 2 , we see that Y = 2 Pi C 

where £ is die damping ratio. The roots of die denominator polynomial are -a±)CO d , 
where CC = £ti) m is the damping coefficient of the m tl ' mode, CP d = 6> m ”\/l - is the 
damped natural frequency of the m th mode, and C0 m = V DA m /p l is the undamped 
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natural frequency of the m ,h mode. Note also that the time constant of the system is 
given by t= 1/a [12]. 

Damping of the beam is achieved, then, directly through the selection of £ or 

indirectly through y. Since Ct= C, C0 m , critical damping occurs when £=1 and a = G> m . 

1/3 

In [19], Maestrello defines a power law relationship a = 0.5 (O d which, unless 
otherwise specified, is the damping coefficient used throughout this study. 

Solving <p k for the m th mode gives 


- 

+ k ~ 


e aA [y sin(co (i A)/p l 


a sin( (ojA) + coj cos(a)jA)]/ coj 


- ( 0 ^ e aA sin(w ( jA)/(Od 


e <xA sin(o»^A)l (O^ 
e~ aA [cos((o c [A) - a sin((OdA)/cod] 


(A .7) 




Appendix B Modal Damping Ratio with Power-Law 
Proportionate Damping 

In [19], Maestrello used a power-law proportional damping of the form 

1. - £>„ * ( B.1) 

where, as discussed in Appendix A, 

(B.2) 

y m is the damping coefficient of the m th mode, C0 m is the natural frequency of the m th 
mode, and is the damping ratio of the m th mode. This is also the form of damping 
used in the computer simulations of Section 9. To determine the effective damping 
ratio, for each mode, Eqs. B.l and B.2 may be combined to form the following 
equation: 


Cl# 



(B.3) 
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Substituting x = CO d yields a cubic equation in X : 


3 X 2 

X + 7 = 0) m , 

4 m 


which has one real and two complex roots. Table 9.1 gives the un 


(B.4) 


damped natural frequency of the first mode as 6)^ = 570.90 rad. /sec. Substituting this 
value into Eq. B.4 and solving for the real root gives X = 68.818. Since 

i = cf'J. ( B.5) 
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r ,,i/3 

yx = co d) 

( B.6) 

and from Eq. B.l, 

= 03 J^~ = 0.007265. 
G), 

( B.7) 

Similarly for mode two, the undamped natural frequency given in 
1573.6 rad. /sec., which gives a real root to Eq. B.4 of X = 135.29. Thus, 

Table 9.1 is COj = 

_ 0.5 yfx 05j^29 _ 0003696 
$2 1573.6 

(B.8) 
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